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Abstract 



Non-magnetic impurities in an S=l/2 Heisenberg antiferromagnetic chain 
are studied using boundary conformal field theory techniques and finite- 
temperature quantum Monte Carlo simulations. We calculate the static struc- 
ture function, S imp {k), measured in neutron scattering and the local suscepti- 
bility, Xi measured in Knight shift experiments. 5i mp (fc) becomes quite large 
near the antiferromagnetic wave- vector, and exhibits much stronger tempera- 
ture dependence than the bulk structure function. \i has a large component 
which alternates and increases as a function of distance from the impurity. 

75.10.Jm,75.20.Hr 



1 



Typeset using REVTeX 



Although the spin chain problem has been a popular topic for theoretical physicists since 
the early days of quantum mechanics, the correlation functions of the antiferromagnetic 
Heisenberg spin- 1/2 chain could only be calculated with the help of modern quantum field 
theory [0]. Adding non-magnetic impurities to a spin-chain compound breaks the chains 
up into finite sections with essentially free boundary conditions. The correlation functions 
in the presence of such a boundary were calculated only rather recently 0. These results 
provide a simple application of a general theory of conformally invariant boundary conditions 
which has been applied to a wide variety of quantum impurity problems in condensed matter 
and particle physics ||. These functions exhibit a universal dependence on the boundary, 
at long distances and times. In this paper we wish to focus on a couple of applications of 
these results of experimental relevance: the impurity contribution to the static structure 
function, S- uap (k, T), and the local susceptibility, Xi(T). We derive analytic expressions for 
these quantities using field theory methods and compare them with finite-T Monte Carlo 
simulations using lengths of up to I = 128 with a varying number of time-steps up to 64 and 
several hundred thousand sweeps through each lattice. 

The Heisenberg Hamiltonian for the antiferromagnetic spin-1/2 chain 

H = JJ2S i -S i+1 (1) 

i 

is equivalent to a free boson field theory in (1+1) dimensions in the low energy, long-distance 
limit The spin operators are expressed in terms of the boson <fi as 

w d x <p/V2^ + a(-l) j cos v 7 ^, (2) 

where a is a constant. The boson Hamiltonian is then simply given by the free part together 
with terms which become irrelevant as the temperature is lowered. Those irrelevant terms 
give rise to temperature and finite length dependent corrections with a characteristic power- 
law. 

This theory has been used successfully to calculate impurity effects 0, the low energy 
spectrum H, and correlation functions Q. The latter agree well with recent neutron- 
scattering experiments ||. Like the expression for the spin operators in equation (0), the 
correlation functions also acquire an alternating and a uniform part as a function of the site 
index x. At finite-temperature, the alternating part is given by: 

< S z {xM)S z {yM) >ait- c^- { ~ ir I (3) 

V(3 ^ ^-y-vM) ginh n( X -y+vAt) 
V V I3 vf3 

(At = t 2 — t\.) We set the lattice spacing to 1. The spin- wave velocity is known to be 
v = Jit/ 2 from the Bethe-ansatz. The constant c can be determined numerically and is 
given by c = a 2 /2 times an arbitrary normalization of the two-point function, which we 
chose to set to one. The irrelevant terms in the Hamiltonian give logarithmic corrections to 
this expression ||. In fact, it has been shown recently that the logarithmic corrections give 
rise to an infinite slope of the uniform susceptibility at zero temperature 

The correlation functions in the presence of a boundary were first calculated in reference 
0. There it was argued that the free boundary condition on the spin operators corresponds, 
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in the continuum limit, to a boundary condition on the bosons: 0(0) = 0l(O) + 4>r(0) = 

y^7r/8. Since 4>l is a function only of vt + x and 4>r of vt — x, this implies that we may simply 
regard this boundary condition as defining cf) R to be the analytic continuation of 0^ to the 
negative axis 4>r(x) = — 0l(— x) + yV/8. Whereas the bulk correlation function factorizes 
into a product of left and right 2-point Green's functions, the boundary correlation function 
becomes a 4-point Green's function for left-movers. Consequently, while the uniform part is 
largely unaffected by an open boundary condition, the alternating part gets modified to 



c (-iy-y- 

v 



sinh^f sinh^ 

vp vp 



(3 \ S inh ^+v+^) S i n h Ax+y-vAt) sinh l^vpM s i n h ^zMzlM 

\ vp vp vp 



(4) 



vf3 



which reduces to equation (Q) in the bulk limit x?/ > \(x — y) 2 — v 2 At 2 \. 

Here we have also included the time- dependence of the Green's function, but we will 
only calculate the equal-time spatial Fourier transform, S(k), deferring consideration of 
the full dynamical structure function to later work. We predict a characteristic impurity 
contribution to the structure factor, which may be observable in magnetic Neutron Scattering 
experiments on quasi one-dimensional spin-1/2 magnetic compounds (e.g. KCUF3). Doping 
with impurities will break the spin-chains and thereby introduce the desired open ends. For 
a finite chain of length / we can define a structure factor Si(k) as 

Si(k) = j < S z {x)S z (y) > e lk{x - y ^ ^ S{k) + . (5) 

x,y=\ 

The structure function for the finite chains has been decomposed into a "bulk" part S(k) 
which is independent of length and an "impurity" part of order 1/7: 

S imp (k) = limZ[S,(*0 -£(*)]. (6) 

The bulk part reproduces the signal of a system without open ends (e.g. an infinite chain) 
while the effect of the open boundary condition is entirely contained in the impurity part. 
Higher order 0(l/l 2 ) terms will also be present, but can be neglected if the impurities are 
dilute. Since each impurity creates the same contribution Si mp (k) in the dilute limit, the 
experimental signal will contain the impurity part as a term which scales with impurity 
concentration n to first order: 

S exp (k) w S{k)+nS imp (k). (7) 

From the results of equations @ and (f|), it is clear that we expect interesting effects for 
wave-vectors near k w tt. Field theory predictions for small k — n and T are obtained by 
Fourier transforming equations (^) and (|j). We assume here that the impurities are dilute 
enough so that the infrared cut-off is always given by the inverse temperature (3 -C l/v . The 
bulk structure function || can then be expressed in terms of the digamma function ip fT0| 
and the reduced variable k' = (k — tt)v/3/ti: 

S{k') = 2c [In (A/3 J) - Re ^ (1/2 - ik' [2)\ , (8) 
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where A is a constant depending on the cut-off. 

The impurity contribution, Sj mp (k') is obtained by Fourier transforming equation (f|) with 
the bulk part, from equation ([3]), subtracted off assuming two open ends. This subtraction 
eliminates the ultra-violet divergence, giving the scaling form: 



Sinip^') = c^- [ dw 

IT JO 



cos k'u I I sinh 2 u \ [°° cos k'u 
du—- — \ 1 r^> 1 — / du- 



o sinh u u sinh w J Jw sinh u 

= v(3f(k>) (9) 

Here u = n(x — y)/vj3, w = tt(x + y)/vj3. Note that, apart from the logarithmic term 
in equation (||), S(k') and Si mp (k') are functions only of the scaling variable k! = {k — 
7r)vf3/7T, but we expect corrections to this scaling behavior from irrelevant operators and 
the ultraviolet cut-off which become smaller as T —>■ and k — it — ► with k' held fixed. 
For small k' we have S(k') oc ln(/5) + const. + 0(k' 2 ) and S imp (k') oc w/3[const. + 0(k' 2 )], 
so that the impurity part has a much stronger temperature dependence. The second term 
in equation <g) can be written / °° dw J™ du = -^limp(l/2 - ik'/2), which is the 

dominant contribution for small k! . At large k', however, the leading behavior comes from 
the first integral Si mp (k') — > —v[3c/irk' 2 , while S(k') vanishes exponentially. (In the opposite 
limit l/v (3 — > oo we get a delta-function at k = it for both the impurity and the bulk 
contributions Si mp (k = ir) oc I and S(k = it) oc ln(/).) 

Our Monte Carlo results for S(k) and Si mp (k) are shown in figures (^]) and (0) for four 
temperatures. Note, that since the impurity part scales with the inverse temperature (3, it 
may make a sizable contribution even at moderate impurity densities. To show the predicted 
scaling form we plotted the results as a function of the reduced variable k' in figures (§) and 
(||). The Monte Carlo simulations agree reasonably well with the field theory predictions. 
This comparison with the Monte Carlo data was used to extract the constant c = 0.14 and 
A = 0.75 in equation (R). 

We now consider the local susceptibility x% a ^ ari Y arbitrary site i under the influence of 
a uniform magnetic field h acting on the complete chain 

= ^ < si > \ h=0 = i E < s;s; >, (io) 

For a chain with periodic boundary conditions, Xi is the same for all sites because of trans- 
lational invariance. 

If we are dealing with an open boundary condition, however, the translational invariance 
is clearly broken and we would naively expect the open end to be more susceptible. Moreover, 
it is now possible, in the field theory treatment, to have a non-zero alternating susceptibility 
as a function of site index Xx = X™ 1 + (""l^Xx*- Using the analytic continuation of the 
left-movers onto the negative half axis from above, Xx* is given by a non-zero three point 
Green's function: 

X f = P<S'^{x) JdySUv)> 

= ^L r dy / ie -i^M*,t>) e iV^M-x,t>)^ { t) + hc \ 

V87T J-oo \ OX 
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dy- 



4tt J-oo s2 s inh ^ + x + iv At) ^ sinh ^(y - x + ivAr) 

(11) 



v M S inh 2™ 

V IT D/3 



where x is the distance from an open boundary condition. At low temperatures the alternat- 
ing part actually increases with the distance yfx from the open end. Any finite temperature 
suppresses this growth exponentially with x, so that we expect a typical maximum which 
gets shifted further into the chain as the temperature is lowered. Furthermore, even at 
T = 0, the staggered magnetization does not increase indefinitely with distance from the 
impurity, but rather oscillates with a wavelength, Airv/h, i.e. 

fo hr 

M alt {x,h,T = 0) = aJ -sm{-£). (12) 

This exotic behavior is similar to Friedel oscillations except that the 1/r 3 decay which occurs 
there gets enhanced to a y/r growth due to a combination of reduced dimensionality and the 
absence of charge fluctuations in this pure spin system. 

The result from equation ([TTJ) can be confirmed independently with quantum Monte 
Carlo simulations. The local susceptibility as a function of distance from the open end is 
shown in figure ([5]) from Monte Carlo simulations at T = J/15. After extracting the uniform 
and alternating parts as shown in figure we can compare the alternating part to the 
predicted form from equation (|TTD where the overall constant was chosen to be a = 0.58. The 
field theory prediction c = a 2 /2, together with the value c = 0.14 from our MC measurement 
of S(k) gives a pa 0.53 in reasonable agreement. While the shape of the theoretical prediction 
for fits the Monte Carlo results very well, there is an unexplained shift of about two 
sites, which might be due to irrelevant operators. The functional dependence in equation (JTTJ) 
holds rather well for all temperatures f3 sampled (up to the shift of two sites). For T = J/15 
the shift in the susceptibility due to the impurity is larger than the bulk susceptibility over 
a distance of about 25 lattice sites from the impurity. Thus we expect that it should be 
possible to observe this effect in nuclear magnetic resonance Knight shift experiments. Note, 
that Xi < for small even i, so that those spins will tend to anti-align with the applied field. 

The uniform part of the susceptibility is not directly affected by the boundary condition, 
but gets an additional non-universal contribution near x = from an irrelevant boundary 
operator which also appears to be present in the Monte Carlo results in figure (0). 

This shift in the uniform susceptibility is what would be expected classically, but the large 
alternating part is a purely quantum mechanical effect. 

In conclusion, we have calculated the effect of impurities on the neutron-scattering cross- 
section and the NMR Knight shift using both field theory and Monte Carlo methods. The 
two methods are in reasonable agreement and the effects seem large enough to be observable 
experimentally. The Knight shift actually increases with distance from the impurity in the 
limit of zero field and temperature. 
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FIGURES 

FIG. 1. The bulk structure factor S(k) according to quantum Monte Carlo simulations. 

FIG. 2. The impurity part of the structure factor Si rnp (k) according to quantum Monte Carlo 
simulations. 

FIG. 3. Monte Carlo results for the shifted bulk structure function, S(k')— 2cln(AJ/3) compared 
to the field theory prediction of equation (||), with c = 0.14 and A = 0.75. 

FIG. 4. Monte Carlo results for the scaled impurity part TSi mp (k') compared to equation @ 
with c = 0.14. 

FIG. 5. The local susceptibility vs. distance from the open end according to Monte Carlo 
simulations at T = J/15. 

FIG. 6. The uniform and alternating parts of the local susceptibility according to Monte Carlo 



simulations at T = J/15 compared to the field theory equation (11) with a = 0.58. 
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FIG. 1. The bulk structure factor S(k) according to quantum Monte Carlo simulations. 




FIG. 2. The impurity part of the structure factor Si mp (k) according to quantum Monte Carlo 
simulations. 
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FIG. 3. Monte Carlo results for the shifted bulk structure function, S(k') — 2cln(A«//3) compared 
to the field theory prediction of equation (8), with c = 0.14 and A = 0.75. 




FIG. 4. Monte Carlo results for the scaled impurity part TSi m p(k r ) compared to equation (9) 
with c = 0.14. 
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FIG. 5. The local susceptibility vs. distance from the open end according to Monte Carlo 
simulations at T = J/ 15. 
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FIG. 6. The uniform and alternating parts of the local susceptibility according to Monte Carlo 
simulations at T = J/15 compared to the field theory equation (11) with a = 0.58. 
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